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We report a study describing the charm quark by a domain-wall fermion (DWF) 
C\| ' in lattice quantum chromodynamics (QCD). Our study uses a quenched gauge 

ensemble with the DBW2 rectangle-improved gauge action at a lattice cutoff of 
a -1 ~ 3 GeV. We calculate masses of heavy-light (charmed) and heavy-heavy (char- 
monium) mesons with spin-parity J p = O^ and IT, leptonic decay constants of the 

o 

f >; ■ charm quark mass is found to be m^ s (m c ) = 1.24(1)(18) GeV. The mass splittings 



X 



charmed pseudoscalar mesons (D and D s ), and the D°-D° mixing parameter. The 



in charmed-meson parity partners A 9) j=o and A g) j =1 are degenerate within statisti- 
Q^' cal errors, in accord with experiment, and they satisfy a relation A q=U( i j > A q=s j, 

■ also consistent with experiment. A C-odd axial vector charmonium state, h c , lies 



22(11) MeV above the Xci meson, or m^ c = 3533(ll) s t a t. MeV using the experi- 



mental Xd mass. However, in this regard, we emphasize significant discrepancies in 
the calculation of hyperfine splittings on the lattice. The leptonic decay constants 
of D and D s mesons are found to be fjj = 232(7) s t a t. (io)chiral(H)syst. MeV and 
/d s //d = l-05(2) sta t. (^2)chirai(2) S yst.) where the first error is statistical, the second 
a systematic due to chiral extrapolation and the third error combination of other 
known systematics. The D°-D° mixing bag parameter, which enters the AC = 2 
transition amplitude, is found to be Be>(2 GeV) = 0.845(24) sta t. (lj 4 ) c hirai(105) sys t. . 
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I. INTRODUCTION 

Properties of hadrons containing one or more charm quarks have been intensively inves- 
tigated in recent experiments . While the main purpose of these 
experiments is to acquire information useful or necessary to CKM phenomenology, it is also 
exciting that several new hadron states have been discovered and confirmed in these experi- 



11 



12 



13]. 



Lattice 



ments. For the latest reviews on these new states, see, for example, Refs. 
QCD should be able to describe various features of these new states. However, in lattice 
calculations of heavy quark systems, the discretization error associated with masses as large 
as the lattice cutoff makes the interpretation of numerical results ambiguous. O ne p ossible 
way to avoid this systematic error is to rely on an effective theory such as HQET |lJ . fisL llfil ] 
or NRQCD 0,^J. With such an approach, however, one must discard either the ability 
to study quarkonia or to take the continuum limit. An alternative is to rely on a relativistic 
method with brute numerical force. Thanks to rapid growth of computational resources, 
discretization errors in this approach for a relatively light heavy quark like charm are begin- 



ning to be 
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jrought under control. Several studies have already been made in this direction 



22, 
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In this work, we explore the feasibility of applying the domain-wall fermion (DWF) 
Hi as a heavy quark within the quenched approximation. Domain-wall fermions have 



cessfully applied to 

H h y , y q 



been successfully appliec 
elements 



to calculations of the light hadron spectrum and weak matrix 
. The primary advantage of working with this fermion action 
is its ability to retain continuum-like chiral symmetry even at a finite lattice spacing, by 
adding a fifth dimension to the lattice. The size of the symmetry breaking is represented by 
the residual mass, am ies , which is adjusted by changing the fifth- dimensional extent of the 
lattice and is of the order of 10~ 3 in lattice units in state-of-the-art calculations. It is known 
that the presence of exact chiral symmetry guarantees the absence of 0(a) discretization 
error, and in the domain-wall case, the violation of the order of am res implies that the leading 
error is significantly suppressed. Thus it seems to be a natural extension to apply the DWF 
formalism to heavier, c and b, quarks. To make the next-to-leading-order discretization 
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error (~ 0((am q ) 2 )) as small as possible, our calculations are carried out on relatively fine 
lattices with a -1 ~ 3 GeV. Although this study is done in the quenched approximation, the 
experience should be useful in the future studies with dynamical fermions. 

One way to test a given heavy quark formalism is to see how consistently the formalism 
describes both heavy-heavy and heavy-light systems with a single value of the mass param- 
eter. While a single mass parameter does not have to describe both systems completely 
consistently in the quenched approximation and at a finite lattice spacing, making a quan- 
titative test at some point will offer a reference for future work. Although we do not take 
a continuum limit here, such a study should become easier in the near future as available 
computational resources grow. 

Once the charm quark mass is determined, it is interesting to see how well the lattice 
calculation reproduces the mass spectrum of recently discovered hadrons: the excited D 
meson states with spin-parity J p = + and 1 + . The spectrum of these mesons had been 
studied before their discoveries by two groups based on chiral quark models These 
turned out to describe some qualitative features well. After the experimental discoveries, 
these groups refined their analyses quantitatively , [3] • The studies suggest that chiral 
symmetry in the light quark sector plays an important role in the spectrum of these mesons. 
Since we apply the domain-wall formalism to both heavy and light quarks, we expect to get 
better descriptions of these mesons than earlier works with Wilson-type light quarks. 

We would also like to make a lattice study of flavor physics, especially weak matrix 
elements of D mesons. As an example, we present the D meson leptonic decay constants 
and the bag parameter, Bd, relevant to D°-D° mixing. Experiments are currently searching 



for evidence of mixing in the neutral D-meson system |5|. Since the expected amplitude is 
very small in the Standard Model, once it has been discovered, it would provide an important 
probe of physics beyond the Standard Model(See, e.g., Ref. 36[.) Unlike the neutral B meson 
system, in D°-D^ mixing it is not clear whether the short-distance contribution mediated 
by a local four-quark operator dominates over long-distance contributions. Even so, it is 
still sensible to have an idea about the size of the short- distance contribution. Knowing 
the size of Bd will be also useful when evaluating Bb by extrapolation in quark mass. In 
this paper we present the first DWF charm lattice study of this quantity. Note there were 



Wilson fermions studies earlier such as Refs. 



0,0. 



The first exploratory study of massive DWF was done in Refs. |22l,|40(, which looked into 
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the low-lying eigenmodes of the five- dimensional Hermitian domain-wall Dirac operator in 
detail and their dependence on the bare quark mass, am q . All eigenmodes are classified into 
one of two kinds of states. The first are the physical, or "decaying" , states which are bound 
to the four-dimensional domain walls located at either end of the fifth dimension. Their 
wavefunctions fall exponentially away from the wall. When am q <C 1, these states dominate 
the low-lying eigenmodes of the whole system and describe four-dimensional physics. The 
other class contains unphysical, or "propagating" , states, which have non-zero momentum in 
the fifth dimension. Their eigenvalues are large, ~ 0(1 /a). From the viewpoint of the four- 
dimensional effective theory, these states are unphysical, a source of non-locality that can 
invalidate the effective theory. The gap between these two types of states is controlled by the 
domain wall height "M5" parameter. In the study of Refs. ^ turned out that as am q 

increases the absolute values of the eigenvalues of the decaying states rapidly increase, and 
the binding to the domain walls becomes less tight or even unbound. On the other hand, the 
eigenvalues of the propagating states increase only slowly. Thus, if am q further increases, at 
some point the lowest eigenmode in the system will be one of the propagating states. Then 
one might worry that something wrong happens to the four-dimensional effective theory. 
In the past quenched studies (39I - I40I ]. this mass threshold can be as low as 0.2 (in lattice 
units) for both j3 = 6.0 Wilson gauge action and (3 = 0.87 DBW2 gauge action, and 0.4 for 
(3 = 1.04 DBW2 gauge action with M 5 = 1.8. In the present work, as will be described, a 
higher lattice coupling is used for the DBW2 gauge action. This results in a threshold of 
am q ~ 0.5 or higher in lattice units, while the bare charm quark mass is below 0.4. 

The rest of the paper is organized as follows. In Sec. [H] we summarize our numerical 
methods and their parameters. Then meson mass spectra and charm mass analyses are 
presented in Sec. IIHI Section |IV] presents results for decay constants. We will also present 
the first DWF calculation of the mixing parameter, Bp, in Section The S'f/(3)-breaking 
ratio is discussed in Section fVTl Extrapolations in l/mh ea v y to the static limit are reported 
in Section IVHI Some systematic errors are discussed in Sec. IVIIII Section |IX] concludes 
with a summary and future outlook for this program. Preliminary results on spectrum and 



decay constants in this work are reported in Refs. 



41 



42 
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II. NUMERICAL METHOD AND PARAMETERS 

The numerical lattice QCD calculations reported in this paper were performed in the 
quenched approximation using the QCDSP computers at RIKEN-BNL Research Center 
and Columbia University. The gauge configurations were generated in a previous RBC work 



kaon mixing bag parameter, Bx |44j . with the DBW2 rectangle 



determining the neutra' 

improved gauge action y, |46( with gauge coupling j3 = 1.22. With the 106 gauge config- 
uration in total used in that study, the lattice cutoff measured from the p meson mass is 
a~ l = 2.914(54) GeV. This implies a physical spatial volume for the 24 3 x 48 lattices of about 
(1.6 fm) 3 . We will use this cutoff estimate as a standard in this report. If we use the static 

n 

quark potential instead, we obtain a^ 1 = 3.07 GeV |44|. Notice that the difference between 
these two determination is about 5%, which is smaller than found on coarser lattices. 

Of the 106 gauge configurations reported in the earlier study, in the present work we use 
only 103, due to the shutdown of the QCDSP computers in early 2006. This results in slight 
differences in the estimations of observables. We confirmed they are not significant. 



We use domain-wall fermions (DWF) |2J, |2Jj, |26j to describe both heavy (charm) and 



light (up, down, strange) quarks. The number of sites in the fifth dimension, L s = 10 and the 
domain-wall height, M 5 = 1.65, are taken to be the same as in the previous work j^J. The 
simulation parameters and some numerical results obtained in Ref. |44j which are relevant 
to this work are summarized in Table |U As seen from the Table, our choice of parameters 
results in a residual mass of O(10~ 4 ) in lattice units, about 0.3 MeV. To study the quark mass 
dependence in a comprehensive way, we use five values each for heavy and light quark masses 
as listed in Table |U They respectively cover approximate ranges of [^m c h ar m, f^charm] and 
[^strange, f^strange] ■ Note that our choice of the heavy mass appears reasonable: Figure [T] 
shows the eigenfunctions in the fifth dimension for the lowest three eigenvalues with 
various masses: 0.03, 0.4, 0.5, 0.7. The largest mass in the present work for the heavy quark 
is 0.5, and as shown in the plot, its lowest three eigenstates drop exponentially away from 
the walls and so do not appear to be unphysical. 

Table HT1 lists the meson operators, ipTip, used in this work. We used source/sink com- 
binations of wall-wall (denoted as C^T 2 ) an d wall-point (denoted as C^ 2 ) types for two- 
and three-point functions. The source position is set at t STC = 7. The gauge fields are 
fixed to Coulomb gauge. The quark propagators are obtained by solving under periodic and 
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anti-periodic boundary conditions in time and averaging over them. This procedure gives 
us a twice larger period in time. In the calculations of two-point functions, we double the 
statistics by using additional set of correlators with t src = 41. 

As was mentioned in the introduction, we calculate the masses of heavy-light (D and D s ) 
and heavy- heavy (charmonium) mesons with spin-parity J p = T and 1 T , decay constants 
of the pseudo-scalar mesons, and the D°-D° mixing parameter. We use the standard single- 
elimination jackknife method for statistical error analysis. Further details of numerical 
methods used in calculating respective quantities are summarized at the beginning of the 
sections reporting the results. 



III. SPECTROSCOPY 



The meson masses are extracted by fitting the wall-point two-point correlators to the 
hyperbolic form 

Acosh(m(t - L t )), (1) 

after shifting the wall source time slice of the correlators to t — 0. 

Using the 103 configurations, we first repeated the measurement of light meson spec- 



troscopy as in Ref. |44j, and confirmed that the differences from that reported in Ref. 44 1 
using full 106 configurations are negligible. For example and for later use, we present the 
bare strange quark mass in lattice units here. We choose to use rriK and m p as inputs to 
determine the strange quark mass. Follo wing the same fit ranges and the same functional 
forms in chiral extrapolation as in Ref. |44J, we obtain am strange = 0.0298(13), which is 
consistent with am strange = 0.0295(14) quoted in Ref. j^jj. 



A. Charmed Mesons 



The effective mass plots for the heavy-light mesons with the four different spin-parity 
states we are discussing are shown in Fig ure 121 for dTfi^^j — 0.4 and o/TJ-Hght — 

0.032, which 

are the closest combination to the bare charm (obtained below) and strange (discussed above) 
mass, respectively. As seen from the figure, the plots for all four mesons show reasonable 
plateaux with this mass combination. Similar quality of plateaux are obtained for the other 
combinations of am heavy = 0.2, 0.3, 0.4 and 0.5 and ami ig ht=0.024, 0.032 and 0.040 as well. 



7 



In addition, for the pseudoscalar (0~) and vector (1~) channels, we could also extract the 
masses for the mass combinations of amheavy = 0.1 and araught = 0.008 and 0.016. We take 
the fit ranges as 15 < t < 23 for 0~ and 1~ states, and 10 < t < 16 for + and 1 + . The 
numerical results for meson masses are summarized in Tables IIII| IIV| fVl and IVII 

We use the pseudoscalar D s meson mass (1968.3(5) MeV), the bare strange quark mass 
estimate of am strange = 0.0298(13) j^J, and p- meson mass m p = 770 MeV as inputs to 
determine the bare charm quark mass. The determination using quarkonium system will be 
discussed later. For this purpose, we first interpolate the heavy-light pseudoscalar masses in 
light quark mass to am stra nge, and then interpolate between the resulting data for amheavy = 
0.3 and 0.4. Through this procedure we obtain an estimate of arnf^^ l p = 0.3583(22), where 
the superscript denotes the inputs used. Using am^™ 1 ', we estimate the masses of the 
four different spin-parity states, resulting in Figure H3 which shows the comparison of our 
estimates (circles) with the experimental values (vertical lines), and in Table IVII II Not 
surprisingly, the calculations are in reasonable agreement with the experiments, to within a 
few %. 

A more stringent test is provided by the parity splittings, A Qt j = m Dq ^j+) — m Dq (j-), 
rather than the masses themselves. The numerical results in physical unit via a" 1 are shown 
in Table IVHI together with their experimental values. First we extract some features from 
the experiments: 

1. A g ~ A g i for both q = ud and s. 

2. A u d t j > A Si j for both J = and 1. 

3. The difference, A udtJ — A S)J , is close to m Da (j-) — mD ud (j-)- This means that the 
positive parity state masses depend on the mass of the light spectator only weakly 
while the negative parity states change by m£, g ^j-) — Tn^jj-) m m s — m u d ~ m s 
between D S (J~) and D u d(J~). 

There have been attempts to understand some of these features using model analyses. As 



we noted in the Introduction, Bardeen et al. 



34j described these charmed heavy-light mesons 



using a chiral quark model in a way that respected heavy quark symmetry, and predicted 
that « A„ with an action that A QCD / mdl „ m corrections are 8m a. U Nowa k 
1. [35] : 



et al. |35( made similar predictions with a slightly different model. Becirevic et al. 



47J noted 
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difficulties in understanding the experimental observation of A u ^ j > A S) j in terms of a 
version of chiral perturbation theory extended to include the four spin-parity states of the 
heavy-light mesons. 

Importantly, many previous lattice calculations have systematically overestimated the 
parity splitting by about 50 to 200 MeV for A sfi 3, H H HQ Q- These works 
are summarized in Figure which is the same plot given in Ref. |53( but the results from 



Dougall et al. 



52j is modified by applying r = 0.50 fm instead of their 0.55 fm. 



Now let us turn to our results. Both A s0 and A s l are overestimated compared with their 
experimental values by 1 - 1.5 standard deviations, which corresponds to a difference of 35 
- 60 MeV. Although our results for A^j are not as precise as A Sj j since the extrapolations 
of the positive parity states to the chiral limit of the light quark suffer from large statistical 
uncertainty (as described below), they show a consistency with the experimental values. 
The degeneracy between A g and A g l are well reproduced for both q = ud and s within 
the statistical uncertainties. A comparison with the previous lattice calculations is made 
in Figure EJ Looking at the data around m c /mg = 1, the present result obtained at our 
relatively fine lattice cutoff turns out to be consistent with that of Dougall et al. [3], which 
is obtained in the continuum limit. 

The heavy quark mass dependence of A s l and A s are shown in Figure El together with 
the available experimental data. The light quark mass is fixed to the strange quark mass. 
It is seen that A Si o and A Sj i are statistically indistinguishable for amheavy > 0.3, but that 
A Si o > A Sj i becomes clearer for smaller omt e!lvy as A Sj o increases while A Si i stays more 
or less constant. The latter behavior may be supported by the experimental values of the 
7^(1270) and K*(892) masses. 

The light quark mass dependence of the splittings is shown in Figure El where am heavy 
is fixed to am chirm ■ Both Aj =0 and Aj =1 moderately increase toward lighter quark mass, 
which is consistent with the experimental observation. 

The splitting between the pseudoscalar D s and D mesons, mn s - mo, is calculated 
to be 101(5) MeV. The central value is slightly above but completely consistent with the 
experimental value of 99 MeV. 

As indicated in Figure El the estimates of the hyperfine splittings are significantly smaller 
than the experimental value. We find the 15* hyperfine splitting vtld* — mz> = 93(4) MeV 
and mo* — mo s = 82(2) MeV, while the experimental results are 142 and 144 MeV re- 
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spectively. This has been a long-standing problem in lattice calculations of heavy quark 

I I 

systems. Take the D s system, for example: UKQCD 54] studied the splitting using a 
nonperturbatively 0(a)-improved Wilson fermion action on quenched Wilson gauge action 
lattices at (3 = 6.2, giving mn* — rriD s = 95(6) MeV using m p to set the scale. In later 
studies, A. Dougall et al.|52j]. using the Sommer scale = 0.55 fm from the K* / K mass 
ratio, obtain rri£>* — m^ s = 97(6) MeV after taking the continuum limit. They also report 
a result of 96(2) MeV from two-flavor lattices at cutoff ~ 1.7 GeV. Di Pierro et al.j^J from 
2+1 staggered dynamical fermion 20 3 x 48, a ~ 0.13 fm (from r = 0.5 fm) lattices report 
the result m D * - m Ds « 112(20) MeV. 

There are many possible reasons for this underestimate in our calculation. Besides the 
quenched approximation, ambiguity in the lattice spacing, and the light quark mass being 
extrapolated from rather heavy values to the physical point, one possible reason is the 
absence of the "clover term". If one requires accuracy through 0(a 2 m charm AQ CD ) for on- 
shell quantities, it turns out that one needs to incorporate the clover term into the DWF 
with proper coefficients j^. Although naively it is expected that the 0(a 2 m c harmAQCD) 
error is fairly small in the present calculation, this needs to be confirmed in future work. 

From the above observations, it is interesting to look at the ratio, 



(2) 



A h yp 

to see how well a given heavy quark formalism describes the whole heavy-light system. Lat- 
tice calculations have overestimated the numerator while underestimating the denominator, 
thus obtaining the ratio closer to the experimental value seems to be extremely difficult for 
any formalism of lattice heavy quark. In our case, we obtain about 5 for amheavy = 0.4, 4 
for 0.3, and 3 for 0.2, and presumably any value of the heavy quark mass would not repro- 
duce the experimental value of 2.4 for the charm-strange system. Although we need to take 
into account the mixing with the other + or 1 + states before a reliable conclusion can be 
drawn, as amheavy becomes smaller, it is worthwhile to keep watching this ratio in future 
calculations. 
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B. Charmonium 

The effective mass plots for the four spin-parity channels of the heavy-heavy meson system 
with amheavy = 0.4 are shown in Figure Similarly good plateaux are observed for other 
heavy quark masses of 0.1, 0.2 and 0.3. We summarize the meson mass estimates obtained 
from them in Table [X] 

If we interpolate to am^-llvm = 0.3583(22) from subsection IIII Al we obtain the charmo- 
nium mass estimates shown in Table IXII and Figure |Hl where the experimental values are 
shown together. Notice that with these inputs the mass of the rj c (J p = 0~) state is consistent 
with the experimental value. Alternatively we can get an estimate of m c harm = 0.3561(11) 
for the bare charm quark mass from the i] c (J p = 0~) and a p . This of course is consistent 
with the am^hl-cm estimate. 

On the other hand the calculated hyperfine splitting of 43(1) MeV is significantly smaller 
than the experimental value of 116 MeV, just like in the charmed meson cases. One of the 
reasons is what we described in the heavy-light case that the lack of "clover" term might 
cause the smallness of the hyperfine splitting. Secondly, heavy quarkonia hyperfine splittings 
are notoriously difficult to reproduce by lattice calculations. In the quenched approximation, 
a detailed study of the relation of hyperfine splitting and the lattice scale was carried out by 
QCD-TARO collaboration [56]. Their result is 77(2) (6) MeV after continuum extrapolation. 
In the dynamical three-flavor staggered case, S. Gottlieb et. al. (stJ reported their hyperfine 
splitting to be 107(3) MeV from the "fine" MILC lattice configurations (a ~ 0.086 fm), 
with lattice scale obtained from bottomonium system. In the charmonium system, the 
problem appears independent of the heavy quark action adopted and is not solved even 
using dynamical configurations. For more details, see the review article in Ref. j^. 

The masses of the P-wave states, in contrast to the charmed meson cases where they 
are overestimated, are several standard deviations underestimated from the experimental 
values. 

Also we note a result of the mass of a yet-to-be-established C-odd axial vector state, 
h c {J pc = e„ta, ly it has been a long - s ta n( W pu ZZ ,e Q. « two 



positive results with a mass of m hc = 3524.4(6)(4) MeV |59l|60| and 3526.2(0. 15)(0.2) MeV 



were reported. A third experiment [62j, though, did not confirm this. In the present 
lattice calculation the mass difference between h c and xa is estimated as 22(11) MeV (see 
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Table IXI|) . This is obtained from the ratio of the h c and x& correlators in the same fitting 
range used for extracting the Xco and Xci masses, so is likely more reliable than the absolute 
values of the excited-state meson mass themselves. Note these masses are significantly 
underestimated in the present work. If we add this difference to the experimentally known 
Xd mass of 3510.59(10) MeV, we get a value of 3533(ll) sta t. MeV. We quote only the purely 
statistical error here. The h c correlator by itself yields 3361(21) stat . MeV, also with only 
purely statistical error. 

Finally we discuss the charm mass and lattice cutoff estimations solely by charmonium. 
First we can determine the bare charm mass by looking at some charmonium mass ratio. 
For example, if we take a traditional choice spin- weighted mass ratios, [{m Xc0 + 3m Xcl ) — 
( m vc + 3 m J/->p)]/ ( m ri c + 3mj^) = 0.13664(6), we obtain an estimation of m c h arm a = 0.235(2). 
Then by matching the experimental mass of (m Vc + 3mj/^)/4 = 3.0676(12) GeV and the 
corresponding interpolation of the calculated mass, 0.804(8) a -1 , we obtain a cutoff estimate 
of a -1 = 3.82(4) GeV. If we use spin-0 mesons alone we obtain m c h ar m« = 0.225(2) and 
ar 1 = 3.90(3) GeV and from spin-1 mesons m c h arm a = 0.239(2) and ar 1 = 3.79(4) GeV. 
This estimate from charmonium is about 30% larger than the one from p meson mass 0]. 
The most likely cause of this is of course the quenched approximation. The rectangular 
improvement of the DBW2 action may also be playing some role here. 



C. Charm Quark Mass 

Presently available estimates of the charm quark mass are given in the PDG review [63] ; 
the charm quark mass is estimated to be in the range 

1.15 GeV < m^(m c ) < 1.35 GeV. (3) 

The renormalized mass can be obtained from 

m ren = Z^ s (m f + m res ) (4) 

where Z^ s = Zj£ is the mass renormalization factor between lattice and continuum. 

Z^ can be obtained with lattice perturbation theory or RI/MOM-scheme non-perturbative 
renormalization either from matching the lattice quark propagator to the continuum one 

or from the scalar bilinear operator renormalization factor as 1/Zg 0, 0. z match 
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matches between MS in the continuum and whatever scheme was used on the lattice; the 
matching between MS and RI/MOM scheme can be obtained from Ref. 0|. Aoki et al.Q] 
calculate one-loop Z^ s with DWF fermions in the DBW2 gauge action, giving Z^ s = 
0.989441; at the scale [i = 2 GeV it gives m c = 1.0330(64) GeV from the D s estimate of 



am 



t charm P= 0-3583(22). Chetvrkin calculates the anomalous dimension of the running of the 
quark mass to 0(af) j3, Q in MS scheme: 

o d 



dfi 2 



-m 



(»/) 



(M) 



m 



(«/) 



(/^)7r. 



("■/) I „ ( n f) 



a* 



(5) 



The renormalization group invariant mass 

"33 - 2N f 



mRi = lim m(/i) 



67T 



(6) 



67 



The running charm mass is then given by 

m c (/i 2 ) = m c , m (a s ) 4/u [l + 0.68733 (a s ) + 1.51211 (a s ) 2 + 4.05787 (a s ) 3 ] 



(7) 



our m^ s (m c 



using a s from Eq. EZJ we find at the scale of m c 
The perturbative renormalization contains a systematic error of 0{a 2 s ) 



1.239(8) GeV. 
4%, giving us 



MS, 



m c ) = 1.24(1)(4) GeV. If we start from the charmonium estimate of m c a = 0.235(2), 
we obtain m^ s (m c ) = 1.07(1) (4) GeV. We find that the difference between these two charm 
quark mass estimates dominates our systematic error; after combining this with the other 
known systematics, we find the charm quark mass to be 1.24(1)(18). 

For the convenience of the reader, some past estimations from quenched lattice QCD 
with continuum extrapolation are listed in Table 11X1 Note, UKQCD recently [73] used a 



two-flavor clover sea to obtain m 



MS, 



rrir 



1.29(7)(13) GeV after continuum extrapolation. 



IV. LEPTONIC DECAY CONSTANTS 



The leptonic decay constants are obtained in a standard procedure from two-point cor- 
relation functions with heavy-light current A4 = \j/;7 M 7 5 \I/ c : 

(0 I At I PS at rest ) = ifps ■ mps . (8) 

We calculated both C^ p and correlators for pseudoscalar mesons. Since these two 
correlators give the same pseudoscalar meson mass (mps), we can extract mps along with 
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amplitudes from fitting these two correlators simultaneously. That is, we minimize the x 2 
given by 

x = 2^ 



t=t„ 



a*4 P (t) 



+ 



C^(t,0)-AZcosh(m PS (t-L t )y 2 



(9) 



where cr(t) is the jackknife error of the correlator at t. The simultaneous fit gives us more 
stable amplitudes for two-point correlators than individual fits for the single correlators. 
The decay constants can be obtained from 

aA a P 

flat = P w _ (IQ) 

where V denotes the spatial volume of the lattice. 

Booth [7^1 and Sharpe and Zhang [3] calculate the chiral behavior of the heavy-light 
decay constants in the quenched approximation to be: 

\fm Qq f Qq = F x + F 2 m q + (F 3 m q + F A ) lnm,, (11) 

after replacement of m 2 oc m q . McNeile and Michael |zf| reported that the effect of the 
chiral log on fg with static-light is small. This indicates that we may be able to drop the 
most divergent term In m q . To be concrete, we enumerate below various fitting functions that 
we will use for the purpose of extrapolating/interpolating the light quark mass dependence: 

Vm Qq f Qq = F x + F 2 m q + F 3 m q \nm q , (12) 
y/m Qq f Qq = Fx + F 2 m q + F 3 m 2 g , (13) 

Vm^fQq = Fx + F 2 m q} (14) 

at fixed amheavy- Then we examine the 1 /amheavy behavior and interpolate to am c harm- 
The upper half of Figure El shows the light quark mass dependence at fixed amheavy of the 
quantity ^mQqfQq extrapolated to — m res according to Eq. Upland Figure [TU] shows the light 
quark mass dependence of the quantity y/mg^fQq interpolated to m stra nge- Then we further 
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linearly interpolate (in 1/M) the heavy quark mass dependence to amh eaV y = aw c h arm , and 
summarize the results in Table IXHI The results from various expressions used for light 
quark extrapolation/interpolation are consistent among the fits. Here we take the fit from 
the quadratic expression, since this falls roughly in the middle of the range from the other 
expressions, and incorporate the difference from the other fits into the systematic error. 
Thus, we obtain 

/jf = 225(7)(^)MeV, (15) 



/g = 243(4)(ljj)MeV, (16) 

before renormalization. 

In order to compare our lattice calculations with continuum results, we need to properly 
renormalize our lattice decay constants by a factor of: 



Z M = Z U x ^DWFKav,) (17) 
y 2 9> DWF (delight) 

The nonperturbative light-light current renormalization in the chiral limit, rrif = — m res , is 
calculated in Ref. Q], as 0.88813(19). Here we employ a quark mass- dependent renormal- 
ization factor 

lat (o) am f (1 + (am f ) 2 ) cosh(mf, 0) ) - 2(am / ) 2 o;sinh(m{ ) 0) ) 

A,,dwf [am f ,u) = (0) , (18) 

(1 — (amf) z ) sinn(mj, ') 

in terms of the tree-level on-shell pole mass 

ifU 2 + ^(1 + am 2 ) 2 + am 2 uo 2 (uo 2 — 4) 



mf ] = In 



-ami 



1 + am j — 2am f co 



(19) 



with uj = 2 — M 5 , and Z^ WF defined as 



QR - \l ^DWF^lat- (20) 



Therefore, the heavy-light current renormalization is 1.0492(22) with m c h ar m — 0.3583(22); 
thus fo = 232(7) (^q) and fo s = 254(4) ( + °) MeV. Recent experimental measurements for 
these values are f D + = 222.6±16.7±|;| MeV 3 and f Ds = 267 ±33 MeV Q and the recent 



lattice three-flavor dynamical results of MILC and Fermilab |78[ are: fjy+ = 201±3±17 MeV, 
and f Ds = 249 ± 3 ± 16 MeV. 
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To calculate the ratio of the decay constants, we first use the chiral form to obtain 

fD s V^~s/fDV^ = 1.11(2)03); (21) 

Figure ITT1 shows an example of our heavy quark mass interpolations of this ratio, with chiral 
interpolation/extrapolation on the light quark mass according to Eq. We may then 
retrieve the ratio of decay constants for the D/D s systems by multiplying in appropriate 
factors of mass: 

W/ D = 1.05(2)(±°). (22) 

The experimental value suggesting the above ratio is 1.19(25) and the MILC number is 
1.24(22). Our result has a smaller central value but agrees with the experimental and 
dynamical result within one a. 



V. D-D MIXING 

In the Standard Model, D°-W mixing is strongly suppressed by CKM and GIM factors; if 
such a mixing is observed, it might be evidence for physics beyond the Standard Model [79] . 
The .D-mixing is the only probe for the dynamics for c-quark. The effective Hamiltonian of 
the QCD contribution to AC = 2 in the Standard Model comes from the box diagram 

H ^ 2 = 0iK:^r (m -~r^ )2 (o + 2C) (23) 

where O = C7 M (1 — 7 5 )wc7 /J (l — j 5 )u and O' = 5(1 + 7 5 )wc(l + 7 5 )u results from non- 
negligible external momentum. The 6-quark contribution is proportional to V u i,mb highly 
suppressed by the very small CKM matrix factor V u b, which is ignored here. The BaBar 
5-factory studies D — D° mixing from the semileptonic decay modes of D* + — ► tt + D° and 
D° -h> Kev [79]. In this work, we focus on the bag parameter of D meson, often defined for 
historical reasons as 

Bd ~ 8^2 f 2 ' ( 24 J 
3'"'DJ D 

where O = C7 M (1 — 75)^07^(1 — 75)^ and Bd is 1 in the vacuum-insertion approximation. 
The parity-even operator for Bp in the continuum limit is 



Ovv+aa = (c7 M M)(c7 Ai w) + (c7 5 7 m m)(c757 m m). 



(25) 
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However, since DWF does not have exact chiral symmetry, Ow+aa mixes with four other 
operators 

Ovv-aa = (c7 m m)(c7 m m) - (c7 5 7 M M)(c757 Ai M), (26) 
Oss±pp = (cm) (cm) ± (c7 5 w)(c7 5 w), (27) 
Ott = {ca^u)(oj^u). (28) 

The mixing coefficients for these operators do not go to zero as mj approaches the chiral 
limit, but are highly suppressed by 0((am res ) 2 ), at least in the case for Bk with DWF. 
This has been theoretically and numerically demonstrated in Ref. Q], and one can simply 
ignore the contribution from the above four operators. However, in the finite quark mass 
region, the correct way of solving this problem is to measure the matrix elements of all the 
operators in Eqns. EBJ EH and EHl directly on the lattice. These operators are relevant to 
beyond the Standard Model sources of mixing in AF = 2 processes |44|. Since we did not 
do this measurement while the data were taken, we quote the maximal uncertainly from the 
mixing coefficients, which are expected to be 0((am^) 2 ), due to the soft chiral symmetry 
breaking from the non-small rrif term. 

On the lattice we can rewrite Eq. Indirectly in terms of the correlation functions obtained 
on the lattice: 



R (lat) <0|x t (^nk)C(t)x t (t Sr c)|0> 
-Dl 



'PS ZnAiP 



Cp w (t, t S nk)Cpw (t, t src ) 



(29) 

snk 



where t src and t sn k are the source and the sink location of the quarks, and the x{t) are quark 
bilinear interpolating meson fields. The results can be found in FigurelT2lwith various mh ea vy 
and miight- The plateau looks pretty good for extracting the value of -Bps? an d details are 
listed in Table IXTTTl 

Booth Q - Sharpe & Zhan g Q c alculate th. chiral behavior of the height eaeson 
bag parameter in the quenched approximation to be: 

Bq q = Bi + B 2 m q + B 3 m q In m q + B± In m q . (30) 

The coefficient of the most divergent term, ln(m g ), is accompanied by a factor of 1 — 3g 2 , 
where g is the coupling of D — D* — n that can be obtained from D* — ► Dti decay 8(|, which 
yields a value of 0.2(7). Note that this g for B — B* — ir is about the same due to heavy-quark 
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symmetry. However, it is still hard to judge the effect of this logarithmic dependence in our 
range of light quark mass. We will ignore it for the rest of paper and incorporate it into the 
systematic error by comparing with the unquenched calculation, where the ln(m q ) term is 
absent. In order to estimate the systematic error due to our choice of fitting form, we adopt 
different expressions to extrapolate/interpolate the light quark dependence: 

B Qq = B t + B 2 m q + B 3 m q \nm q ; (31) 



B Qq = B 1 + B 2 m q + B 3 m 2 q , (32) 

B Qq = B 1 + B 2 m q . (33) 

Here we show the extrapolations using the Eq. HH1 in the upper graphs in Figure to 
the chiral limit, mj — — m res and in Figure El to the strange quark mass. The results from 
various fits are summarized in Table IXIVI and different analyses agree with each other within 
statistical error bars. Therefore, we take the quadratic fit as the central value and absorb 
the other fits into the systematic error; thus we have 

B% = 0.859 (24) (tf) (34) 



Sg = 0.848(7)(ig) (35) 

before proper renormalization. Note that since D s is a charged meson, there cannot be 
oscillations between D s and D s . Here we calculate its bag parameter merely for the interest 
of studying the SU(3) breaking effect in charmed systems. Also, the extrapolation to the bag 
parameters from B and B s mesons are simply side products of this study. It is interesting 
to check the goodness of 1/M extrapolation with the static quark lattice studies. 

We take advantage of the bag parameter nonperturbative renormalization (NPR) done 
in Ref Q with RI/MOM scheme and further convert into MS. Ref. HQ calculated the 
conversion factors in NLO perturbation theory. First, we convert the renormalization factor 
into an RI value 



J- H T— -'RI/MOM 



z la t ,RI/MOM ( ^ lat) (36) 
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where 

a s (/i) 



12tt 



i 918 - 90N f - 2AN] lnln(/x 2 /A| CD ] 



(33 - 3N f ) ln(/i 2 /A 2 QCD ) ^ (33 - 2iV» 2 ln(/i 2 /A| CD 

with N f = 0, and A° CD = 238 MeV and Jri/mom(^/ = 0) = 2.883; that gives Zf 
1.409(5). Then we calculate the factor for Nf = 4, as 



(37) 



'MS 



25/3 



1 + T- 

4:71 



MS 



(38) 



The Aqcd with N f = 5 is 217 MeV |6j, suggesting that the Aq C d is 276 MeV for Nf = 4. 
Then a s (fi = 2 GeV) is 0.283. Therefore, Z^ s , is 1.000(4) and our bag parameters are: 



B D = 0.845(24) (± 24 ) 



(39) 



J B Ds = 0.835(7)(+ 2 ), (40) 

ZW^ = 0.987(22)(1° 7 ). (41) 

VI. SU(3) BREAKING RATIOS 

In the long run, lattice QCD should provide a high-precision determination of SU(3) 
flavor-breaking of the AF = 2 heavy-light matrix element 

M Qs = (PqsIQj^I - 75)3^7^(1 - 7 5 )g|PQ 8 ) {42) 
Mqi (P Ql \Q lfl (l- l5 )lQ lfl (l- l5 )l\P Ql ) ' 

where Q stands for heavy quark, I stands for light quark (u,d), and Pqi represents a pseu- 
doscalar meson composed of Q and I. There are two ways of obtaining this ratio: first, by 
calculating the matrix element directly [3]; secondly, by calculating the decay constants and 
bag parameters separately and combining them with Eq. In this work, we will focus on 
the second method (often referred to as the "indirect" approach) to get the ratios by the 
combination of ratios of bag parameters and decay constants: 

_ M Qs _ B Qs (m Qs fQ S ) 2 
Q M.qi B Ql (m Ql f Qi y 1 ' 
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Therefore, it is important to obtain the ratio of 

¥\f¥=^ < 44 > 

JQi y a Q l 

Table IXVI summarizes £q and the SU(3) breaking for different chiral extrapolation formulae. 
From the charm quark sector, we obtain 

£ c = 1.071(23X1°!), (45) 

and 

r c = 1.273(65)(ljj 7 ). (46) 



VII. EXTRAPOLATION TO STATIC B MESONS 



Although our main goal in this work is to apply the DWF to the charm quark physics 
directly, it is of some interest to extrapolate the heavy quark mass to the static limit. 



In the past, as is pointed out in Ref. 



84L l85| . some static limit B parameters obtained 



from extrapolation from charm region with simple linear function in l/m^eavy were found 
disagreeing with direct static calculations: after all the charm mass may not be sufficiently 
heavy to justify such a simple extrapolation to the static limit. However, as a mere by- 
product of our work on charm, the extrapolation may still be instructive. 



A. Decay constants 

Since we do not know how to renormalize the heavy-light decay constant for a static 
quark with DWF as the light quark action, we will only focus on ratios of decay constants. 
Since the bottom quark is so heavy compared to the scale of our physics, we may extrapolate 
to the static quark limit to get the result pertaining to B mesons. The ratio is 

tf B .y/^J fay/mSf*** = 1-08(4) (47) 

Taking m^/m^ as 1.017 we have 

(/^//B) static = 1.06(4)(^). (48) 

A previous quenched study using the static approximation for the heavy fermion action and 
the step-scaling technique |86j gives 1.11(16) when extrapolated to the continuum limit, 
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which is consistent with our extrapolation result here. A recent s 



tudy using partially 

quenched two-flavor DWF lattices with the static approximation [87[ (with lattice cutoff 
~ 1.7 GeV) yields 1.29(4) (4) (2) for the same quantity, giving us an idea of the systematic 
error due to the quenched approximation. HPQCD use a 2+1 staggered fermion sea with 
NRQCD heavy quark action and lattice cutoffs w 1.6 GeV and ~ 2.6 GeV to obtain the 
ratio 1.20(2)(1) |8j|. 

B. Bag parameters 

Following the analyses in Sec. El we extrapolate lattice bag parameter to the static limit 

(j3*f) static = 0.940(33) (if ) (49) 

^tjstatic = o.919(9) (ij|). (50) 

We set the scale /x at the mass of tVquark, 4.5 GeV, and the renormalization factor Z|f s (m&) 
(with Nf = 5) is 0.920(3), which suggests 

(B B ) static = 0.865(33) (if) (51) 



(^jstatic = . 845 (9)(+3) (52) 



Ref. 



89j uses static heavy and Wilson light quark actions to get B B (rrib) = 0.98(4) (if 8 ) 
at « 1.8 GeV lattice cutoff. Another quenched study with static heavy and overlap light 
quarks gives 9(3| (again, with the finest lattice cutoff ~ 1.8 GeV) 

(B Bs (m b )) static = 0.940(16) (22). (53) 

Our number does not agree too well with previous static results. This might be due to the 
coarse lattices used in their simulations. However, our number agrees better with JLQCD's 
quenched calculation with NRQCD 0, where the finest lattice spacing is 2.3 GeV: 

B Bd {m b ) = 0.84(3)(5), (54) 
B Bs /B Bd = 1.020(21)(^)(j|), (55) 
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B s .(m 6 )=0.85(l)(5)(iS), (56) 

after taking the continuum limit. Our results for the B case are consistent with the two-flavor 
DWF static-light calculation in Ref. Q: 

(jBi?) static = a812(48)(67)( +O oo)) (57) 

{B B .)«« 3 ° = 0.864(28) (71) (±° 20 ), (58) 

(B Bs /B B )^ tic = 1.06(6)(3)(1), (59) 

and the two-flavored 0(a)-improved Wilson fermion sea with NRQCD heavy quark study ji^ 
done by JLQCD 

(B B ) = 0.836(27)(j£)(B J ,./fl s ) = 1.017(16) (tf 7 ). (60) 

C. SU(3) breaking ratio 

Extending the discussions in Sec. IVI| we have the results for the bottom sector: £& = 
1.019(37) (±34) and the SU(3) breaking ratio r b = 1.274(94) (j^,). It is useful to recall that in 
the first calculation in Ref. [83J, using a Wilson fermion action, found = 1.30(4) (til) and 
1.17(2)(^g 2 ) after linear continuum extrapolations from direct and indirect methods, respec- 



87]. 



tively. A previous static quark study on a two-flavor DWF sea gives: £& = 1.33(8) (8) 
Using JLQCD ratio of bag parameters J3] and HPQCD 2+1 result [z^l on decay constants, 
the Lattice '05 review [9J| gives & = 1.210(1^) and the previous world average 94] eives 
^ = 1.23(6). Our central value is smaller than these previous results but quite consistent 
with other studies given our relatively large errors. We first take the light quark mass limit 
according to the heavy-light meson bag parameter chiral formula and then take the heavy 
quark mass to rn^fc^] Figures fTTTl and ITU show that the chiral behavior of the light quark 
part makes the bag parameters for mesons with up and down quarks somewhat larger than 
those with a strange quark. However, the other studies get their bag parameters by fitting 
the bag parameter as a function of heavy-light pseudoscalar mass, instead of the heavy-light 
chiral forms; see the dependence for our data in Figure El The bag parameter increases 
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with the pseudoscalar mass. Therefore, if one fits the bag parameter function as input of 
heavy- light pseudoscalar mass, from a relation ttl^b,d} s / m {B,D} > 1, which always holds, it 
follows that -B{b,d} s /-B{b,d} > 1- Therefore, £q in these studies must come out larger than 
our values. 



VIII. SYSTEMATIC ERRORS 



In the previous sections, we present our calculation with statistical errors which are 
computed using the jackknife procedure, along with systematic errors mostly caused by 
chiral extrapolation. There are other potential sources of systematic errors caused by the 
quenched approximation, finite-volume effects, operator mixing and matching calculations. 
Here is a detailed estimate of the various systematic errors. 

Finite volume: The proper way to estimate this error is to perform the calculation on 
at least two different volumes and extrapolate to the infinite- volume limit. However, 
in this work, we only perform our calculation in one lattice volume, which is about 
(1.6 fm) 3 . Since the scale of D system is in between K and B, we will estimate our 
finite-volume error by quoting whichever system (K or B) has larger finite-volume 



error in previously published studies. Ref. 



44j . which uses the same lattice box we do, 



quoted a 2% error in B^ by comparing with a 2.4 fm box. We also compare our Bb 
found by extrapolation to the static limit with a quenched, larger- volume (~ 2.4 fm), 
calculation j^jj, which gives BB d (rrib) = 0.84(3) (5) with NRQCD fermions; this gives 
us an estimation of 0.4% error. Taking the larger of these two estimates, we conclude 
that the finite- volume effect on the Bp should be smaller than 2%. 

Continuum extrapolation: In this paper, we only perform our calculation at one lattice 
spacing, and therefore we cannot extrapolate our result to the continuum limit. This 
should be checked carefully in future works. The next-best thing we can do is to 
estimate the error from a continuum extrapolation study from the same gauge config- 
uration. Such a study has been performed on the decay constants and matrix elements 
of the K system in Ref. j^J with an additional lattice cutoff at 1.982(30) GeV. The 
result shows mild dependence from a -1 m 3 GeV to the continuum. Therefore, we add 
0.2% to our systematic errors due to continuum extrapolation. 
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RI formulation: We need the a s to first convert the scale-dependence RI/MOM NPR 
factor, 2^ I//MOM to an RI value and further convert it to an MS one. In the expression, 
we use the NLO order formulation, which leaves the remaining leading error up to 
0(a 2 s ) « 4%. 



NPR renormalization factor in bag parameters: We adopt the RI/MOM NPRrenor- 
malization factor from Ref. [44], and we quote the estimation within that paper as 1%. 

Operators mixing: We have discussed mixing with wrong chirality operators, which con- 
tributes around 10% uncertainties to our final B D calculation. 

Quenched approximation: The quenched approximation ignores sea quark loop contri- 
butions which in general are considered to be a major contribution to systematic 
errors. Ref. |4J|, which uses the same quenched lattice ensemble as in the present 
work, quotes a 6% error on the Bk factor due to this approximation, after comparing 
with the number calculated on two-flavor DWF lattices. Similarly, we can com par e 
our static quark limit value with the one calculated on two-flavor DWF lattices (87 1, 
which gives us 6%. Therefore, it seems reasonable to use a 6% systematic quenching 
error for Bp as well. 



IX. CONCLUSION 



In this work, we use the domain-wall fermions (DWF) formulation for charm quark as 
well as the three lighter flavors, up, down and strange, on the lattice with a relatively high 
cutoff (around 3 GeV). We use the mass ratio m^/m p to set the bare strange quark mass 
^strange^ = 0.0298(13) in lattice units. Then combining this and another hadronic mass 
ratio mD s /m p we obtain the bare charm mass, m c ^ arra a = 0.3583(22). 

Using these bare quark mass values we found, we conclude the following for charmed and 
charmed-strange meson states: 

• the masses of the J p = ± and 1 ± D, D*, D s and D sJ states are well reproduced to 
within a few percent. 

• Their parity splitting, Aj, are better reproduced than previous works, with only 10- 
20% over estimations. 
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• The experimental observation of A ud > A s is reproduced. 

• The hyperfme splittings are only 60-65% reproduced. 
Regarding the dependence on heavy quark mass, 

• Aj =0 and Aj =1 are degenerate for mh eavy a > 0.2 — 0.3. 

• Aj =0 increases as mhcavy decreases further, while 

• Aj =1 does not. 

Also with the bare charm quark mass we worked out the renormalization factor, Z m , from the 
existing one-loop domain- wall perturbation calculation. After RG scaling, the m^ s (m c ) = 
1.239(8) GeV. The perturbative renormalization contains a systematic error of 0{a 2 s ) ~ 4%, 
giving us m*®(m c ) = 1.24(1) (18) GeV. 

In the charmonium system we find the i] c mass agrees well with the experimental value. 
The J ftp mass (hence the hyperfme splitting) is smaller than the experimental one, con- 
firming the long puzzle in lattice QCD. We also note a prediction for the mass of yet- 
to-be-discovered C-odd h c state: the mass difference between h c and xa is estimated as 
22(11) MeV. This would translate to the mass of 3533(ll) stat .MeV for this meson. 

The leptonic decay constants are also calculated. We tried to estimate the systematic 
uncertainty from the associated quenched logarithm by adopting three different fitting for- 
mulations. Our result after considering other systematic uncertainties due to the quenched 
approximation and continuum extrapolation give us 



We also discussed extrapolation to the static limit in terms of l/mh ea vy which gives results 
consistent with previous static calculations. 

The bag parameters in the D and D s (purely theoretical but interesting in regard of SU (3) 
breaking effect) are studied. The use of domain-wall fermions on the current fine lattice with 
only softly broken chiral symmetry gives us some advantages such as absence of complicated 
mixing and the availability of the RI/MOM nonperturbative renormalization techniques. 



/D = 232(7)(lg)(ll) MeV 



(61) 



fn./fn = 1-05(2)(1°)(2). 



(62) 
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We include a detailed estimation of systematic uncertainties. The biggest systematic errors 
come from the quenched approximation and the mixing of wrong-chirality operators. Our 
result is: 

B D {2 GeV) = 0.845(24) (^ 4 ) (105), (63) 

iW#D = 0.987(22) (1° 7 ) (23), (64) 

where the first error is statistical, the second systematic from fitting, and the third combining 
all other known systematics. Thus a minor S'f/(3)-breaking is seen in our calculation. We 
also discussed extrapolation to the static limit in terms of l/m hcavy which suggests the D 
and D s meson results are reliable. 

In conclusion, using DWF to simulate the charm quark on a quenched ensemble at a 
moderately high cutoff of about 3 GeV obtained with the DBW2 action results in reasonable 
descriptions for most of the calculated meson observables: the masses and their splittings, 
leptonic decay constant, and AC = 2 mixing. It seems the charm quark propagation is 
successfully described with the current set up. This obviously is an attractive direction to 
proceed, especially with dynamical QCD ensembles that are becoming available. 
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TABLE II: Meson states in this study are created by local operators of the form ipTip in this table. 
Their spin and parity are also listed, as well as charge conjugation for quarkonium cases. 



r 2S+l L 7 j pc 



75 


'So 


o-+ 


li 


3 s 1 


1— 


1 


3 Po 


0++ 


757i 




1++ 


7;7j 




1+- 
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TABLE III: Heavy-light J p = meson mass obtained from fitting to the hyperbolic cosine form 
in the range of 15 < t < 23. 

"Iheavy "Might ™PS ^ Per d.O.f. 



0.008 0.3254(15 

0.016 0.3373(13 

0.1 0.024 0.3495(12 

0.032 0.3617(11 

0.040 0.3739(11 



0.008 0.4715(17 

0.016 0.4812(14 

0.2 0.024 0.4913(13 

0.032 0.5016(12 

0.040 0.5120(11 



0.008 0.5916(19 

0.016 0.6003(15 

0.3 0.024 0.6095(13 

0.032 0.6189(12 

0.040 0.6285(11 



0.008 0.6917(21 

0.016 0.7000(16 

0.4 0.024 0.7086(14 

0.032 0.7175(13 

0.040 0.7266(12 



0.008 0.7734(22 

0.016 0.7815(17 

0.5 0.024 0.7898(15 

0.032 0.7984(14 

0.040 0.8072(13 



0.22 
0.15 
0.09 
0.08 
0.07 



0.36 
0.30 
0.23 
0.19 
0.18 



0.51 
0.46 
0.37 
0.32 
0.31 



0.71 
0.70 
0.59 
0.52 
0.49 



0.84 
0.92 
0.84 
0.78 
0.75 
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TABLE IV: Heavy-light J p = 1 meson mass obtained from fitting to the hyperbolic cosine form 
in the range of 15 < t < 23. 

"iheavy "Might ™PS X 2 per d.o.f. 



0.008 0.4083(45 

0.016 0.4153(35 

0.1 0.024 0.4223(30 

0.032 0.4310(26 

0.040 0.4401(24 



0.008 0.5259(37 

0.016 0.5324(28 

0.2 0.024 0.5400(24 

0.032 0.5484(21 

0.040 0.5572(19 



0.008 0.6304(34 

0.016 0.6373(25 

0.3 0.024 0.6450(21 

0.032 0.6532(19 

0.040 0.6618(17 



0.008 0.7203(32 

0.016 0.7275(24 

0.4 0.024 0.7351(20 

0.032 0.7433(18 

0.040 0.7518(16 



0.008 0.7943(31 

0.016 0.8017(23 

0.5 0.024 0.8095(20 

0.032 0.8176(18 

0.040 0.8260(16 



0.55 
0.40 
0.27 
0.20 
0.17 



0.65 
0.52 
0.36 
0.27 
0.23 



0.63 
0.56 
0.41 
0.33 
0.29 



0.59 
0.57 
0.44 
0.36 
0.33 



0.64 
0.68 
0.57 
0.50 
0.47 



33 



TABLE V: Heavy-light J p = + meson mass obtained from fitting to the hyperbolic cosine form 
in the range of 10 < t < 16. 

™heavy flight ™PS X* Per d.O.f. 

0.024 0.656(16) 0.01 

0.2 0.032 0.654(17) 0.01 

0.040 0.658(10) 0.01 

0.024 0.762(13) 0.02 

0.3 0.032 0.760(10) 0.02 

0.040 0.764(8) 0.02 

0.024 0.854(14) 0.01 

0.4 0.032 0.853(11) 0.01 

0.040 0.857(9) 0.01 

0.024 0.934(14) 0.01 

0.5 0.032 0.932(10) 0.01 

0.040 0.936(8) 0.01 
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TABLE VI: Heavy-light J p = 1 + meson mass obtained from fitting to the hyperbolic cosine form 
in the range of 10 < t < 16. 

™heavy flight "IPS ^ Per d.O.f. 

0.024 0.689(17) 0.004 

0.2 0.032 0.688(12) 0.002 

0.040 0.692(10) 0.002 

0.024 0.792(16) 0.011 

0.3 0.032 0.790(12) 0.006 

0.040 0.793(10) 0.005 

0.024 0.881(16) 0.017 

0.4 0.032 0.879(12) 0.010 

0.040 0.881(9) 0.008 

0.024 0.955(16) 0.02 

0.5 0.032 0.952(12) 0.01 

0.040 0.955(9) 0.01 



TABLE VII: Summary of mass splitting results (in MeV). 

Experiment This work 

A 444(36) 533(90) 

Ai 420(36) 452(78) 
1- - 0- 140.64(10) 93(4) 

A 50 348.4(9) 411(40) 

A S i 345.9(1.2) 380(37) 
1- - 0- 143.8(4) 82(2) 
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TABLE VIII: Summary of the heavy- light spectrum (in MeV). 



Meson (J J 


Experiment 


This work 


D (0 ) 


1869.4(5) 


1867(4) 




2010.0(5) 


1961(4) 


T~\$= / r\-\~ \ 

£>o(° ) 


2352(50) 


2401(89) 


D' l {l+) 


2427(26)(25) 2413(76) 


Df(0-) 


1968.3(5) 






2112.1(7) 


2051(2) 




2317.4(9) 


2379(40) 




2459.3(1.3) 


2431(37) 



TABLE IX: A list of the past charm quark mass estimations from quenched lattice QCD with 
continuum extrapolation. 

Group m^ IS (m c ) GeV action 

Kronfeld [95] 1.33(8) Clover 

Hornbostelet al. [96] 1.20(4) (11) (2) NRQCD 

Becirevic et al. [9_fJ 1.26(4) (12) Clover 

Juge [98] 1.27(5) Clover 

Rolf and Sint [99] 1.301(34) Clover 

de Divitiis et al. [19] 1.319(28) NRQCD 

Nobes et al. [1QQJ 1.22(9) Fermilab 

This work 1.24(1)(18) DWF 
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TABLE X: Summary of quorkonium spectrum in lattice unit. 

(""heavy 0"+ 1~ 0++ 1++ 1+" 

0.1 0.4614(12) 0.5113(24) 0.604(11) 0.632(9) 0.656(17) 

0.2 0.7118(10) 0.7392(14) 0.824(7) 0.848(6) 0.858(8) 

0.3 0.9226(8) 0.9406(11) 1.026(6) 1.046(5) 1.053(7) 

0.4 1.1008(8) 1.1129(9) 1.200(6) 1.218(6) 1.227(6) 

0.5 1.2467(7) 1.2541(8) 1.350(6) 1.364(7) 1.377(5) 



TABLE XI: Summary of the charmonium spectrum (in MeV). The last row presents a separate 
estimate on the mass difference between h c and Xci obtained directly from themesons propagator 
ratio: it is likely more reliable than the calculated excited mesons mass themselves which are 
underestimated. 

Charmonium (J PC ) Experiment This work 



Vc(0- + ) 

JMi-) 

Xco(0 ++ ) 

Xcl(l ++ ) 

h c (l+-) 



2980(1) 
3096.916(11) 
3415.19(34) 
3510.59(12) 



2987(12) 
3030(11) 
3282(21) 
3336(21) 



not established 3360(21) 



h c - Xd 



22(11) 



TABLE XII: Fit expressions and various decay constant related values, 

fit a^y/mEh a^y/m^fp. ^ V^IbT^ (a^^fgj^ (^g^) 8t 

a + bm q 0.0621(18) 0.0677(13) 1.090(19) 0.0718(31) 0.0773(25) 1.073(56) 

a + bm q +cm 2 q 0.0621(18) 0.0688(13) 1.107(19) 0.0718(31) 0.0783(21) 1.085(37) 

a + bm q + cm q In m q 0.0630(21) 0.0688(13) 1.092(26) 0.0725(39) 0.0783(21) 1.076(53) 
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TABLE XIII: List of the B$ values, 
a^hcavy amiight am PS 

0.2 0.016 0.792(10) 0.4839(14) 

0.024 0.792(7) 0.4940(12) 

0.032 0.795(6) 0.5042(11) 

0.040 0.799(5) 0.5146(11) 

0.3 0.016 0.830(11) 0.6033(15) 

0.024 0.828(8) 0.6124(13) 

0.032 0.830(6) 0.6218(11) 

0.040 0.834(6) 0.6314(11) 

0.4 0.016 0.855(12) 0.7027(15) 

0.024 0.853(9) 0.7113(13) 

0.032 0.855(7) 0.7202(12) 

0.040 0.858(6) 0.7293(11) 

0.5 0.016 0.874(13) 0.7865(23) 

0.024 0.871(9) 0.7943(20) 

0.032 0.873(7) 0.8024(18) 

0.040 0.876(6) 0.8110(17) 



TABLE XIV: Fit expressions and and various bag parameter values. 
fit Bp B Db BpJBp (B B ) static (£? Bs ) static (B Bs /B B ) stl 

a + bm q 0.845(16) 0.849(8) 1.001(12) 0.923(22) 0.921(10) 0.998(14) 

a + bm q + cm\ 0.859(24) 0.848(7) 0.987(22) 0.940(33) 0.919(9) 0.977(28) 
a + bm q + cm q In m q 0.874(33) 0.848(7) 0.970(32) 0.958(45) 0.919(9) 0.959(39) 



TABLE XV: Fit expressions and various bag parameter values. 

fit /ds / Bp s r fp s m Ds \ 2 Bp B f Bs j B Bs / /fl s m.B s \2 B Bs 
fp V B D v m D f D > B D [b V B B V m B f B > B B 

a + bm q 1.064(29) 1.257(73) 1.019(53) 1.27(13) 

a + bm q + cm 2 q 1.071(23) 1.273(65) 1.019(37) 1.274(94) 
a + bm q + cm q In m q 1.048(31) 1.219(80) 1.001(53) 1.23(13) 
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FIG. 1: Fifth-dimensional s-dependence of the three eigenvectors of Dh with smallest eigen- 
values (displayed in right- and left-projection pairs): \^ 2 S \ (with y-axis in log scale) with m/ € 
{0.03,0.4,0.5,0.7}. The propagating states only arise for arrif > 0.5. 
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FIG. 2: Effective mass plots with "close-to-quark-mass" simulation points: amneavy = 0.4 and 
amiight = 0.032. 
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FIG. 3: The spectrum of the D s (above) and D (below) systems. The circles are our results with 
statistical errorbars and the horizontal lines correspond to experimental data with one o error. 
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FIG. 4: The summary plots of the previous lattice estimate for A s j including ours. 
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FIG. 7: Effective mass plots in the heavy-heavy sector at the simulation point am-Heavy = 0.4 
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FIG. 8: The spectrum of the charmonium system. The circles are our results with statistical 
errorbars and the horizontal lines correspond to experimental values. 
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FIG. 9: The light (above) and heavy (below) quark mass dependence of a 3 / 2 s/mofD with the chiral 
expression: ^rnQ q fQ q = F\ + Fim q + {F^m q )\n.m q and linear (in 1/mneavy) fit respectively. The 
top figure displays fixed amjieavy G {0.2, 0.3, 0.4, 0.5} from top to bottom. In the bottom figure, the 
triangle (blue) point represents the physical D meson point, 0.0630(21), and the pentagon (red) 
point represents the static quark limit point of 0.0725(39). 
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FIG. 10: The light (above) and heavy (below) quark mass dependence of a 3 ^ 2 ^/mu s fo s with the 
chiral expression: yJrriQ q fQ q = F\ + Fifriq + {F^m q ) lnm g and linear (in 1/rriHeavy) fit respectively. 
In the bottom figure, the triangle (blue) point represents the physical D s meson point, 0.0688(13), 
and the pentagon (red) point represents the static quark limit point of 0.0783(21). 
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FIG. 11: The heavy quark mass dependence of the SU(3) breaking ratio 



/m D f D 



- with chiral 



interpolation/extrapolation f Eg. 112(1 in light quark mass (black points). The triangle (blue) point 
represents the physical charm quark point, 1.092(26), and the pentagon (red) point represents the 
static quark limit point of 1.076(53). 
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FIG. 12: The time dependence of pseudoscalar meson bag parameters at fixed light quark mass: 
0.016, 0.024, 0.032, 0.040 from top to bottom. In each subgraph, we display the heavy quark mass 
dependence of -Bps- 
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FIG. 13: The light (top) and heavy (bottom) quark mass dependence of -Bps with the chiral: 
Bq q = B\ + B^rriq + B^m q lnm q and linear (in 1/M) fit respectively. The top figure displays light 
quark mass interpolation to — m res at fixed on-heavy 1 0.2, 0.3, 0.4, 0.5 from top to bottom. In the 
bottom figure, the triangle (blue) point represents the physical D meson point and the pentagon 
(red) point represents the static quark limit point. 
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FIG. 14: The light (top) and heavy (bottom) quark mass dependence of -Bps with the chiral: 
Bq q = B\ + B^rriq + B^m q lnm q and linear (in 1/M) fit respectively. The top figure displays light 
quark mass interpolation to m s t ra ge at fixed afn-heavy : 0.2, 0.3, 0.4, 0.5 from top to bottom. In the 
bottom figure, the triangle (blue) point represents the physical D s meson point and the pentagon 
(red) point represents the static quark limit point. 
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FIG. 15: The B% as a function of the pseudoscalar meson. The circle points are the data points 
and the black line is the fit of the form: -Bps(m-ps) = B\ + B^w^^ + B^rripg In rap S . The star (blue) 
point is the B l jj point obtained from Fig |13l the triangle (red) point is the B^f from quadratic fit, 
and the pentagon (purple) point is the B 1 ^ from linear fit. 



